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Abstract 



The LHCb RICHl detector uses hybrid photon detectors (HPDs) as its optical sensors. A calibration system has been constructed 
to provide corrections for distortions that are primarily due to external magnetic fields. We describe here the system design, 
construction, operation and performance. 
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1. Introduction 

The LHCb experiment 1 1] is dedicated to precision measure- 
ments of CP violation and rare decays of B hadrons at the Large 
Hadron Collider at CERN. Key measurements of the LHCb 
physics program are outlined in [2J. Excellent discrimination 
among different particle species is vital to study decays with 
similar topologies and to suppress combinatorial backgrounds. 

LHCb has two RICH detectors to provide hadron identifi- 
cation in the momentum range 2-100 GeV/c. Both RICH de- 
tectors use arrays of hybrid photon detectors (HPDs) to detect 
the Cherenkov photons created in their radiators |[Tl[3]|4l. The 
photon detectors need to provide accurate photon position mea- 
surements. This is a critical component of the Cherenkov angle 
resolution, which determines the separation power of different 
particle species, particularly at high momentum. Indeed, at all 
momenta, good photon position resolution is required by the 
likelihood function that is typically used in particle identifica- 
tion algorithms. 

The HPD used in the LHCb RICH detectors |5 6| is an op- 
toelectronic imaging vacuum tube that has a quartz window of 
spherical section, an S20 multi-alkali photocathode, and a sil- 
icon pixel sensor. Its internal electrodes are operated at 16-20 
kV, setting up an electric field cage in which the photoelectrons 
are cross-focused, and follow long (~125 mm) drift trajecto- 
ries from the photocathode to the pixel sensor The pixel sensor 
is bump-bonded to a custom binary readout chip. There are 
32x32 effective pixels in the sensor, each of area 0.5x0.5 mm^. 
However, each effective "LHCb" pixel actually consists of eight 
physical "ALICE" pixels of area 0.5x0.0625 mm^. The eight 
ALICE pixels are grouped together in the readout to form a sin- 
gle LHCb pixel (termed "pixel" herein). The nominal electro- 
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Static demagnification factor of five gives a measurement gran- 
ularity of 2.5x2.5 mm^ at the exterior of the quartz window. 

HPDs are well known to be sensitive to external magnetic 
fields Q [S |9l. In particular, they are efficient only below 
1.5 mT axial fields, and less so if the field is off-axis. The mag- 
netic Lorentz force deflects photoelectrons from their nominal 
trajectories causing pinwheel and shift distortions of the entire 
image and hence a subsequent loss of accuracy in the detected 
photon position. Depending on the angle of the magnetic field, 
this distortion can be quite severe (up to several pixels), com- 
pletely compromising the photon position resolution. 

Both LHCb RICH detectors are affected by the fringe mag- 
netic field of the spectrometer magnet, but it is most pronounced 
in the RICHl detector which is located immediately upstream 
of the magnet. The field is reduced by an iron shield surround- 
ing RICHl (seen in Fig. [TJ, and further by cylindrical Ni-Fe 
alloy shields around each HPD. Even so, initial estimates indi- 
cated up to 3.0 mT fields at the HPDs. With possible inhomo- 
geneities and saturation effects, the resulting non-uniform fields 
could produce serious distortions which would be impossible to 
predict in an a priori manner The calibration system described 
in this paper is dedicated to mapping and correcting these mag- 
netic distortions in RICHl I.10J . A different system is employed 
for this purpose in RICH2 ifTTIl . In addition to magnetic distor- 
tions, there are electrostatic distortions possible from variations 
in the applied HPD high voltages, and optical distortions from 
the quartz window. Distortions from all sources are combined, 
and this calibration system will correct for these effects as well. 

In the next section, the overall design of the magnetic dis- 
tortion calibration system (termed MDCS herein) is provided, 
followed by a description of the relevant details of the system, 
its testing and installation. Subsequently, the data-taking proce- 
dure and analysis method are developed, then the main perfor- 
mance metrics and results are presented. Concluding remarks 
follow. 
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Figure 1: Schematic view of the RICHl detector. The MDCS system sits in the 
photon funnel, which is located just outside of the gas enclosure. There are two 
MDCS systems, one for each of the upper and lower HPD aiTays. Each HPD 
array has 7 columns of 14 HPDs each. 

2. System Design 

The basic requirement on the MDCS system is to remove 
the potentially severe distortions from the array of HPDs. The 
corrected photon position should contribute an error small in 
comparison with the intrinsic spatial resolution of the HPD. 

The strategy for the MDCS is implemented in two steps. 
First, a light spot is projected at a precisely-known position on 
the quartz window of a given HPD in the array. The spot has a 
small size in comparison to the pixel size, and is oriented front- 
on to minimize refractive effects at the window and shadowing 
by the magnetic shields. The spot is moved to scan the HPD 
array in two dimensions, with small enough step sizes to map 
out the distorted image of the pixels. The data taken from this 
scan form a "direct" mapping of the distortion at a given mag- 
netic field. Second, an "inverse" mapping is created, assigning 
a given distorted pixel hit to its real position on the photocath- 
ode. This may be done by a parameterized functional fit to the 
data if the distortion is sufficiently smooth to do so with a small 
enough reconstruction error, or alternately by a look-up table 
if the distortion is too severe or non-uniform. With proper sur- 
vey information, this strategy can provide an absolute positional 
calibration of each pixel in the HPD array. 

Implementing this strategy for the MDCS system required 
consideration of the tight physical constraints on the system in 
terms of overall area (1283x540 mm^), form factor (<40 mm 
high), and maintenance of a clear photon aperture for the HPD 
array. Additional considerations were made due to operation in 
a magnetic and radiation environment. 

The basic system design is shown schematically in Fig. |2] 
Instead of a single light spot moving in two dimensions, spatial 
constraints on the system size forced a solution having a long 
bar with many individual light sources (LEDs), termed a "light 
bar." Each LED in the light bar can be individually powered. 




Figure 2: MDCS overall system drawing. In this view, HPDs are located above 
the system, facing down (into page). Their positions are denoted by dotted cir- 
cles. The light bar movement is horizontal, and it is shown in its fully-retracted 
(parked) position. 

and each is collimated in order to project a small light spot onto 
the HPD focal plane. The light bar is mounted as a gantry be- 
tween two linear-motion control stages. The movement of a 
given light spot in the direction of the stages is very fine and 
nearly continuous, while the positioning of the light spot in the 
direction along the light bar is discrete, effected by powering 
adjacent LEDs in turn. This choice is dictated by the spatial 
constraints on the system. 

3. System Description 

Two identical and independent MDCS systems were de- 
signed and constructed for RICHl, one for each of the upper 
and lower HPD enclosures. In this section, the description per- 
tains to only one of the two systems. 

3.1. Light Bar 

The light bar consists of two long carbon fiber side rails on 
which are mounted a series of 19 LED PCB units, custom de- 
signed for this applicatiorj^ 

Each LED unit consists of a six-layer PCB, on which are 
mounted four LED arrays, an on-board addressable microcon- 
troller, a driver circuit to control each LED, and a collimator 
unit. This unit is shown in Fig. [3] 

The LED array^each consist of a 5 x7 matrix of green LEDs 
on 2.54 mm centers. The LEDs have dominant wavelength 569 
nm and half-width 30 nm. The four LED arrays are tilted as a 
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Figure 3: A close-up of the light bar design, showing a full LED PCB unit, as 
viewed from the HPD array. 

unit at an angle of sin"' (0.5/2.54) = 11.4°, in order to provide 
a smaller effective separation of light sources. This tilt gives 
a 0.5 mm granularity of LED centers, as projected in the long 
dimension of the light bar. During construction, the LED arrays 
were precisely aligned to reference holes in the PCB, epoxied in 
place, then soldered. The PCB was populated afterwards with 
the remaining components. 

The microcontrolle:[^on each PCB is addressed individually 
and downloaded with a pattern code via an SPI bu^which runs 
along the light bar To power an LED, a microcontroller output 
line goes high, sourcing the current for one or more selected 
rows of five LEDs. Concurrently, the microcontroller fires one 
or more column power mosfet£] to sink the current. The 
pattern code selects which row(s) and column(s) are activated. 
All necessary LED array patterns may be activated with this 
logic. Each LED PCB may be programmed with a different 
pattern or all with the same pattern. 

The collimator unit is an array of individual collimators, one 
for each LED in the LED array, shown in Fig. |4] Each colli- 
mator is centered on an LED, and consists of a length of open 
tube machined in a block of Delrin® acetyl homopolymer. The 
collimator geometry has two 6.35 mm lengths of 0.343 mm di- 
ameter separated by one 12.7 mm length of 1.0 mm diameter. 
Extensive design studies of the collimator geometry were car- 
ried out, with the projected spot size as the figure of merit. The 
spot size was designed to be 6.5 mm FWHM at the effective 
HPD focal plane, 85 mm from the collimator 

A practical consideration was the ability to reliably and pre- 
cisely machine an array of 2590 identical collimator holes in a 
light bar. This influenced the decision to use Delrin®, which is 
suited for CNC machining in its ability to be cut without bind- 
ing the tools, and which has low moisture retention hence good 
dimensional stability over time. The collimator unit was con- 
structed as a stack-up of two pieces, due to the limitations in 
the thickness of material though which a small-bore hole can 
be drilled rehably. These two pieces were visually inspected 
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Figure 4: An exploded view of a single LED PCB unit. Light is emitted from 
LED array L into the collimator unit, encountering first part B having short 
small-diameter holes followed by longer larger-diameter holes, and then part A 
having the same short small-diameter holes. Light emerges upwards toward the 
HPD array (not shown). 

after manufacture, and re-worked by hand as necessary to as- 
sure each small hole was free of any residual obstructions. 

Mechanically, the collimator unit is the precision reference 
point for fabricating the light bar. The collimator unit was as- 
sembled on the populated LED PCB using the same reference 
holes as the LED arrays. A precision jig was used to align all 
19 completed LED PCB ^ in the light bar by locating the exit 
faces flat on a reference plane and transversely by using these 
same PCB reference holes transferred to the face of the colli- 
mator units. The result was a ~1200 mm long, planar, two- 
dimensional array of 2590 collimator holes precisely located 
on 0.5-mm projected centers, across LED and PCB boundaries, 
with reference marks for surveying into RICHl. 

3.2. Motion Control 

Motion control was effected by two linear stepper motor 
stage^ with a long enough travel to cover the entire active area 
of the HPDs. The stepper motors allow for precision move- 
ment of the light bar An optical position encoder gives feed- 
back position information for additional accuracy. Both stages 
are wired to a single motor controller, synchronizing the move- 
ment of both stages. Each stage was modified with a 90-degree 
helical gear box to fit the tight spatial constraints. The gearbox 
also allows for manual turning from outside of RICHl, should 
the stages fail. 

Each stage is fitted with two limit switches which indicate the 
home and end-of-travel positions. During movement, a swing 
mechanism keeps the light bar cabling from binding. Mount- 
ing pins allow the removal and replacement of the stages with 
precision. Similarly, the light bar is mounted with locating pins 
from one end for precision reference. 

3.3. Master Control Board 

The light bar and motion stages are both controlled by the 
master control board, a custom PCB which has a master micro- 
controller, internet interface, LVDS drivers, power drivers, and 



The last LED PCB in the light bar is half-sized, containing only two LED 
arrays. All other 18 LED PCBs are full-sized. 
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Figure 5 : A photograph of one MDCS system, installed in a photon funnel prior 
to installation in RICHl. The light bar is seen extended at ~ 1/3 of its full travel. 
During RICH operation, it is folded away into its parked position. 

AC relay. The master microcontrolle:[^was programmed with 
its own small operating system to respond to specific commands 
arriving from the internet porj*^ This allows all commands to 
be issued via TCP/IP. Depending on the command, the master 
control board will either pass the command to the LED PCBs, 
pass it to the stepper motor controller, or act on the command 
itself. Hence by synchronizing the LED pattern and movement 
of the light bar, a cartesian grid pattern (or other required pat- 
tern) may be projected on the entire HPD array. 

The light bar commands are transmitted to the LED PCBs 
over a long cable using LVDS logic to reduce noise pickup, 
and received by a transition board at the head of the light bar. 
Here the signals are level-adapted to the TTL SPI bus which 
communicates to the slave microcontroller on each LED PCB 
individually or collectively. 

The master control board can also turn off the DC power to 
the light bar and AC power to the motor controller. This func- 
tion is quite important, since even when idle but powered, these 
devices may be a source of noise for the HPDs during operation. 

4. Testing and Installation 

Both MDCS systems were tested extensively for mechanical 
and optical performance. 

The motion control stages have a position resolution of 50 
yum bidirectionally. However, a repeatability of ~10 fim is ob- 
tained when traveling unidirectionally (i.e., when the movement 
includes no backlash in the lead screw). So in operation, any 
such effect is reduced by always approaching a new position 
for the light bar from the same direction. 

The Ught output of each spot was set via a current-limiting 
resistor to have a 1 MHz count rate as measured by the HPD 
at the center of the spot, which is equivalent to the highest rate 
expected in actual operation. 



'"Microchip microcontroller PIC 24FJ128GA006. This is the master to the 
slave microcontrollers which are on the LED PCBs. 

"Lantronix XP1002001-03R XPort SE embedded ethernet server 



Figure 6: A pixel hit map of a single LED light spot (with ALICE pixels in the 
y direction). The black dot marks the center of a 2D Gaussian fit, and the cross 
marks the position of the arithmetic mean, for comparison. The smattering of 
other hits are from background processes, as discussed in the text. 

Both MDCS systems were mounted in their respective pho- 
ton funnels before installation in RICHl, as shown in Fig. |5j 
The positions of the light bar at several points along its travel 
were surveyed with respect to the photon funnel. After installa- 
tion, the photon funnel was surveyed with respect to the RICHl, 
and thence defined with respect to the LHCb spectrometer. 

5. Operation in RICHl 

Special data-taking runs (called MDCS scans) are taken pe- 
riodically for monitoring purposes, or in the event a new cali- 
bration is needed (e.g. HPD replacement). Data are taken with 
both polarities of the spectrometer magnet (at 1 .5 T nominal) as 
well as with the magnet off. 

The light bars are controlled interactively via TCP/IP using a 
PVSS projecp] integrated into the RICHl global control sys- 
tem. The LEDs are synchronized with light bar movement to 
create a grid pattern of light spots over the entire HPD array. A 
scan is divided into "steps," uniquely characterized by a fixed 
bar position and a given configuration of LEDs. For simpUc- 
ity, all LED PCB units are configured identically, with only one 
LED per unit switched on. 

During a scan, the light bar moves 42 times and illuminates 
each HPD column, in turn, for ~5 sec. Typically, the sparse 
LED pattern results in an HPD being illuminated by a single 
LED at each step. In this way, the pixel hits can be directly 
associated with a specific LED position. For a given light spot, 
several pixels are hit. The hit pixels are grouped together into a 
cluster, called herein a "peak" (as illustrated in Fig.|6]|. 

6. Calibration Data Analysis 

The first tasks in analyzing calibration scan data are to iden- 
tify the peaks (clusters), and to clean the images in order to 
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Figure 7: The radial width a of cleaned peaks, for magnet on data (solid red) X (LHCb pixels) 

and magnet off data (dashed blue). 



distinguish peaks corresponding to LED light spots from those 
due to unwanted background processes. 

6.7. Peak Finding 

The peak-finding algorithm forms peak candidates from the 
photoelectron hits occurring during the 6000 triggers that are 
recorded at each step. Pixels with more than 100 hits/step and 
with one or more of the eight neighboring pixels having no hits, 
are flagged as noisy, and removed from the analysis. 

Peak candidates are formed in a square window of 7x7 LHCb 
pixels where the central pixel contains at least 1 % of the total 
hits in the HPD (for that particular step) and the neighboring 
pixels have fewer entries. Overlap between these matrices from 
different peaks is not allowed. Selected peak candidates are 
then fitted using a two-dimensional Gaussian function. 

The peak width distributions of the width of the peaks for 
magnet off and on data are compared in Fig. [7] and show that the 
peak shapes are not significantly distorted by the magnet. The 
average peak radius is 0.84 pixels which corresponds to a light 
spot of 5.5 mm FWHM at the photocathode, to be compared to 
the expected design value of 6.5 mm FWHM. 

6.2. Sources of Background 

Three independent sources of background contribute fake 
peaks in the scan data: ion feedback [TT], total internal refrac- 
tion within the photocathode window, or noisy pixels in the sil- 
icon chip. 

Ion feedback occurs when a photoelectron ionizes residual 
gas atoms in the HPD tube. The ion drifts to the photocathode 
producing secondary electrons which are focused on the sen- 
sors with a typical delay of 200-300 ns. Usually the fake peaks 
due to ion feedback are located near the chip center This back- 
ground occurs only in certain tubes in the HPD array. 

Backgrounds from total internal reflection occurs whenever 
an LED illuminates the periphery of the quartz window. Light 
entering the tube may be reflected back into the quartz window 
by the flat surface of the HPD dynode. It is then incident on 
the external surface of the window at a sufficiently large an- 
gle for total internal refraction to occur The curvature of the 
window dictates the radii at which the reflected-then-refracted 
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Figure 8: Fitted peak patterns in an HPD image, (a) Typical image, before 
cleaning, with both signal and background peaks visible, (b) Same image, after 
cleaning, and before distortion correction. Only signal peaks remain. 



light forms ancillary peaks. This background occurs at a similar 
radius, regardless of azimuthal angle, and produces the circu- 



lar features visible in Fig. 8(a) This background occurs in all 
HPDs. 

HPDs with noisy pixels or those affected by ion feedback 
can be easily identified since they always have hits, even when 
not illuminated. These tubes usually have a total number of re- 
constructed peaks that is well above the average of other tubes, 
integrated over the entire scan. 

6.3. Pattern Cleaning 

A pattern-cleaning algorithm is used to distinguish primary 
signal peaks from background arising from the above sources. 
Four variables are used to characterize the properties of the 
peaks and to build a likelihood that optimally selects signal. 
These variables are: (1) the height of the fitted 2D peak, (2) 
the number of hits recorded across the HPD in which the peak 
is found, (3) the proportion of these hits that contribute to the 
peak, and (4) the displacement of the peak in the direction of 
the light-bar's travel, compared to the average displacement of 
all the peaks found in the other HPDs in the column. 
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An unbinned maximum likelihood fit is performed on these 
four variable distributions and combined into a likelihood for 
the signal and the background hypotheses. The logarithmic 
difference between the two likelihoods (Delta Log Likelihood: 
DLL) is taken as the discriminant variable for the peak selec- 
tion. The DLL distributions for single and non-single peaks is 
shown in Fig. |9] As expected, HPDs which have a single fit- 
ted peak have a high likelihood of being a primary signal peak. 
The cut in DLL can be varied depending on the purity required, 
however usually a DLL > cut is applied. If more than one 
peak passes the DLL cut, the one with the highest DLL is kept. 
With this technique, a large fraction of useable peaks are re- 
tained from HPD images containing multiple peaks by requir- 
ing DLL > (HPD images with a single peak are unambigu- 
ous and always retained). An example of the cleaned pattern is 
shown in Fig. |8(b)| The resulting set of cleaned peaks is used 
to extract the final calibration parameters. 

7. Calibration Extraction 

The distortion of the photoelectron image in the presence 
of the magnetic field is observed to be modest and smoothly 
varying. This allows an inverse mapping strategy, mentioned 
above, to be implemented via a parameterization of the distor- 
tion rather than a look-up table procedure. Parameters are ob- 
tained for each HPD and stored in an LHCb database for use in 
the Cherenkov ring reconstruction with collision data. 

7.7. Parameterization 

Due to design of the RICH and its shielding, the magnetic 
field inside an the HPDs is, to a good approximation, parallel 
to the tube axis. The absence of a significant transverse field 
simplifies the parameterization, as radial symmetry about the 
tube axis can be assumed. The transformation of the position 
of a peak on the silicon chip to that of the corresponding light 
spot at the exterior of the HPD window involves both radial and 
axial components. Any point on the pixel chip has coordinates 
{Ra,^a), viz., the radius and azimuthal angle with respect to 
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Figure 10: MDCS distortion correction applied to a typical HPD: (a) before 
distortion correction, and (b) after distortion correction. The cleaned peak po- 
sitions at the photocathode before and after correction are shown as dots. The 
expected grid patterns from two LED units ai'e superimposed. 



the center of radial symmetry on the chip (chip center). Simi- 
larly, the corresponding point on the outer surface of the quartz 
window has coordinates (Rc, (pc)- The transformation is param- 
eterized as 
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Figure 9: DLL distribution of the peak fitting, showing single peaks (blue cir- 
cles) and non-single peaks (red triangles). 



Rc = piRA+PiRl+PiRl 

<Pc = (PA + eo + eiRA+O2Rl+03Rl 



(1) 



where p, and 0, are the parameters to be extracted. The radial 
part is a combination of the electrostatic magnification, the op- 
tical refraction at the curved surfaces of the photocathode win- 
dow, and the magnetic distortion correction. The axial part de- 
scribes the characteristic "S-shaped" displacement, as is evi- 
dent in Fig. |8] These parameters are expected to depend only 
on the magnetic field strength and to be stable in time. The 
whole photocathode image (and so, the chip center) can move if 
small amounts of charge build up on insulating surfaces within 
the tube. This is measured periodically by fitting a circle to 
the uniformly-illuminated HPD disk image from proton-proton 
collision data. This procedure automatically takes into account 
misalignment of the chip with respect to the tube axis. 
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8. Performance on Calibration Data 
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8.1. Position Resolution 

The distribution of residuals between the distortion-corrected 



peaks and the expected LED position is shown in Fig. 12 Clear 
benefit is seen from application of the fitted parameterization, 
particularly in the ALICE pixel direction, which has the higher 
granularity. The variance of the post-transformation data can 
be interpreted as the residual uncertainty on the position of in- 
cident photons due to the magnetic distortions. A simple Gaus- 
sian fit gives a width cr x 0.49 mm (outside the HPD quartz 
window) which corresponds to an factor of ~3.7 improvement 
with respect to the uncorrected distribution. This remaining un- 
certainty compares well with that due to the binary resolution 
of the LHCb pixels (0.81 mm). So it can be concluded that the 
residual magnetic distortions inflate the effective pixel contri- 
bution to the Cherenkov angle resolution by just ~17%. 



Figure 1 1 : Polynomial fits of distortion pairs AR and A(/> as a function of Rc- 



7.2. Fitting Procedure 

The calibration parameters are deduced by minimizing the 
distance between the position of the selected peaks, suitably 
transformed with Eq. [T| and the cartesian grid that is defined 
by the geometry of the light bar and the logic of its movement. 
Once cleaned, only one peak per HPD remains; so it is trivial 
to associate it with the originating light spot, and hence provide 
strong constraints in the fit. As the LED unit is shorter than 
the diameter of the HPD tube, an HPD is typically illuminated 
by two LED PCB units which result in two distinguishable grid 
patterns. Depending on the step number within the LED illumi- 
nation sequence, the light may come from one unit or the other, 
but not from both in the same step. Therefore these two distinct 
grids are easily separated in the data, as seen in Fig. [TO] The 
spacing between light spots at the photocathode and the spacing 
between two grids (inter-LED unit distance) is imposed using 
the design specifications of the light bar. The absolute positions 
of the light spot grids in the xy plane is not imposed but fit- 
ted together with the distortion parameters to allow for a global 
misalignment of the HPD tube. 

Each peak position and corresponding grid point form a 
"pair". The pair is described by the parameters ARc and A^c^ 
orthogonal distances between the two points when projected 
onto a 2D surface perpendicular to the cylindrical axis of the 
HPD. Pairs with ARc and/or AcpcRc greater than twice the spac- 
ing between two consecutive LEDs are not considered in the fit, 
providing additional protection from residual backgrounds. An 
example of the relation of ARc and A(pc versus radius Rc is 
shown in Fig. [TT| These shapes are parameterized by simple 
polynomial functions f7"8l, and then combined with a coeffi- 
cient of magnification to arrive at the full transformation, Eq.[T] 
The peak positions after the full transformation are shown in 
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Fig. 10(b) and exemplify the validity of the parameterization. 
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Figure 12: Typical position resolution of corrected peak positions in calibration 
data, as measured at the HPD quartz window, for (a) x (LHCb pixel) direction, 
and (b) y (ALICE pixel) direction. The uncorrected peak positions are indicated 
by the dashed histogram. 



8.2. Magnification 

The demagnification power of the HPDs is extracted from 
these data and compared with the electrostatic design expec- 
tation of a factor five. In the absence of magnetic distortions 
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9. Performance on Collision Data 
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The most important performance metric of the MDCS system 
is its effect on the RJCHl Cherenkov angle resolution in proton- 
proton collision data. For saturated (fi - I) Cherenkov rings in 



collision events with hadrons in the final state. Fig. 14 shows the 
per-photon Cherenkov angle resolution for the same data, with 
and without the MDCS coiTections applied. With the MDCS 
corrections, a fit with a single Gaussian plus a second order 
polynomial background gives a resolution of 1.597 ± 0.002 
mrad. This is close to the expected value of 1.55 mrad from 
simulation. Without the MDCS corrections, the resolution is 
about a factor of two worse, as determined by a comparison of 
the FWHM. 
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Figure 13: Magnetic distortion correction to the magnification for different 
scans and HPD liigh voltages over a two year period, for a typical HPD. (a) 
Radial correction, (b) Axial correction. 



10. Conclusions 

The LHCb RICHl magnetic distortion calibration system 
described here is dedicated to correcting the photon positions 
from the effects of magnetic fields and other distortions in the 
HPD photons detectors. The MDCS system has been in oper- 
ation since the start of LHCb data-taking. An analysis method 
has been developed to extract distortion correction coefficients 
from calibration data. Applying these corrections to proton- 
proton collision data, the per-photon Cherenkov resolution an- 
gle resolution is improved by a factor of two, and now ap- 
proaches the expectation from simulation. This improvement 
greatly enhances the particle identification performance of the 
LHCb experiment, and indicates that the MDCS system has ful- 
filled its design requirements. 



(using data taken whilst the spectrometer magnet is off) the ef- 
fective demagnification of the HPDs, including the fish-eye re- 
fraction of the curved photocathode window, is measured. 

A second order polynomial is sufficient to parameterize the 
magnification plus refraction effects. The total magnification 
factor, averaged over all HPDs, is 5.6. The linear term domi- 
nates, having an average value 5.3, with the quadratic term be- 
coming significant at high Rc, near the perimeter of the quartz 
window. 

8.3. Correction Stability 

The stability of the MDCS corrections over time is related 
to HPD behavior and long-term MDCS performance. Fig. [13] 
shows the axial and radial magnetic distortion corrections (with 
magnification constant) for a particular HPD at different times 
and high voltage conditions over nearly a two year period. For 
the entire Rc range, the radial and axial corrections are stable 
and within the resolution of the method. This legitimates the 
use of the same calibration parameters during several months of 
experimental operation. A new calibration is performed in the 
event of HPD replacement or significant changes in operation 
(e.g., change of high voltage). 
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Figure 14: The per-photon Cherenkov angle resolution in 201 1 collision data: 
(a) before MDCS distoilion correction applied, and (b) after MDCS distortion 
correction applied. 
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